A New Type of Traveling Interface Modulations in a Catalytic Surface Reaction 
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Pattern formation in reaction-diffusion systems covers 
a wide range of fascinating phenomena in liquid phase 
chemistry, biochemistry, biology and catalytic surfaces 
[1-3]. In general, the patterns arise due to the coupling 
of a non-linear reaction term with diffusion. Reaction 
fronts, target patterns and spiral waves, stationary con- 
centration patterns and chemical turbulence have been 
seen. Various additional factors like global coupling, dif- 
fusional anisotropy, energetic interactions and cross diffu- 
sion of reactants may add to the complexity and diversity 
of the chemical wave patterns. 

Extended bistable systems generically exhibit fronts 
(also called interfaces or domain walls) connecting one 
phase in one part of the spatial domain to the other 
phase in some other part of the domain. In two spa- 
tial dimensions the most natural geometry is a straight 
line for the front position, suitably defined as some inter- 
mediate level curve of the solution. However, already in 
simple bistable systems, initially straight interfaces be- 
tween two domains may undergo a number of instabili- 
ties, see, e.g., [4, Chapter 2] for an overview. A typical 
case is a linear transverse instability leading to a regu- 
lar (periodic) or irregular bending of the front, but with 
small amplitude, which may then often be described by 
Kuramoto-Sivashinksky type of equations, see [5]. An- 
other possibility is that an instability does not saturate 
at some small amplitude, which may yield "fingering" 
and labyrinthine patterns [6-8]. Similar wave instabili- 
ties also occur in excitable media, see, e.g., [9]. 

Here we report on a new type of instability and self- 
organization of an interface, namely interface modula- 
tions that originate from corners and travel along the 
interface in a pulse like fashion, leaving the interface posi- 
tion almost unperturbed behind. These excitations have 
been observed in the NH3 + O2 reaction on a Rh(llO) 
surface. The effect is attributed to diffusional mixing of 



two spatially separated adsorbates causing an excitability 
which is strictly localized to the vicinity of the interface 
of the adsorbate domains. Combining a bistable with an 
excitable system, we set up a general model which re- 
produces the traveling interface modulations seen in the 
experiment. 

The reaction we study is the catalytic ammonia ox- 
idation with O2 on a Rh(llO) surface under low pres- 
sure conditions (10~ 5 mbar) in a UHV chamber equipped 
with a photoemission electron microscope (PEEM) as 
spatially resolving method. Illuminated with a D 2 dis- 
charge lamp (5.5-6 eV) photoelectrons are ejected which 
allow an imaging of the local work function with a spa- 
tial resolution of w 1/zm and the temporal resolution of 
video images (20 ms). At elevated temperatures (T>400 
K) both reactants dissociate upon adsorption into their 
atomic constituents O a d, N ad , and H a(i [10, 11]. The 
atomic adsorbates recombine, forming N 2 , NO, and H 2 
as main products. Also, H 2 is produced and desorbs at 
a high rate, and hence the coverage Oh is always small. 
The adsorbates N and O form a large number of ordered 
reconstruction phases on Rh(110) but under our reac- 
tion conditions only the (2xl)-N/(3xl)-N correspond- 
ing to 9n = 0.5/0.33 and the c(2x6)-0 corresponding to 
0q = 0.66 are relevant [12]. In addition, a mixed coad- 
sorbate phase c(2x4)-20,N may form. 

Over a broad range of parameters the reaction exhibits 
simple bistability, i.e. one observes a broad hysteresis in 
the reaction rates in heating/cooling cycles. The reactive 
branch is associated with the (2xl)/(3xl) of nitrogen, 
the unreactive branch with the c(2x6) of oxygen. Tran- 
sitions between the two states occur via fronts. If one 
adjusts conditions close to equistability both phases are 
simultaneously present as shown by the PEEM image in 
Fig. la. Since oxygen adsorption strongly increases the 
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FIG. 1: Experimental observation of interface excitations in 
the NH3 + O2 reaction on Rh(llO). Experimental conditions: 
T=740 K, p(NH 3 )=3.85 x l(T 5 mbar, p(0 2 )=1.35 x 1(T 5 
mbar. (a) PEEM image showing the interface between oxy- 
gen covered (dark) and nitrogen covered surface area (bright). 
The inset representing an enlarged view of the interface re- 
gion near S shows the formation of dark boundary layer at 
the interface within the oxygen phase, (b) Enlarged view of 
the region marked in (a) showing the pulse-like propagation 
of an interface modulation, (c) Position vs. time plots show- 
ing the temporal variation of the interface position. The data 
were taken from cross sections perpendicular to the interface 
at points E and S in (a). 

work function (WF) (A$ max <~ 1.0 eV) high O a d cover- 
ages are imaged dark whereas adsorbed nitrogen which 
only causes a maximum WF increase of 280 meV appears 
bright [13]. 

The position of the interface is nearly stationary but 
one notices small lateral displacements of the interface 
which emanate near the sharp corner in the phase bound- 
ary and then propagate in a pulse-like manner along the 
interface. This process is depicted in more detail by the 
frames in Fig. lb displaying an enlarged section of the 
PEEM image in (a) . The velocity of the pulse- like exci- 
tations is about 6 /im/s. Cross sections of the interface 
showing the temporal variation of the interface positions 
at two different points, E and S, arc displayed in Fig.lc. 
Near the sharp corner, in point S, the amplitude is below 
the detection limit. Further away, at point E. the am- 
plitude is substantial varying between a few /im and 20 
/im. One notes a drift of the average interface position 
of about 15 /im over an observation time of 170 s which 
is due the fact that the equistability conditions arc not 



exactly met. The time series exhibits irregular behavior 
but the excitability of the interface is quite stable and 
can be observed over several hours. The average period 
of the local excitations is around 10 s. In our experiments 
we found no correlation between the interface angles and 
interface excitations, and the crystallographic directions 
of the surface. 

In order to understand why the excitations remain lo- 
calized at the interface and do not extend into the in- 
terior of the phase it is helpful to look into the chemi- 
cally rather similar system Rh(110)/NO + H2 which can 
be considered as well understood [13, 14]. Some spec- 
tacular chemical wave patterns including rectangularly 
shaped target patterns and spiral waves and traveling 
wave fragments were found there. The excitable behav- 
ior in this system was shown to be based on a cyclic 
change of three different structures; the c(2x6)-0 of oxy- 
gen, the (2xl)/(3xl)-N of nitrogen and the c(2x4)-20,N 
as mixed coadsorbate phase. In the NH3 + O2 reaction 
only two of these three structures are present as stable 
phases while the mixed coadsorbate phase is missing. Ap- 
parently the mixed phase does not form by coadsorption. 

If we assume that by surface diffusion this mixed phase 
may form, its formation is restricted to a boundary layer 
along the interface where the two separated adsorbates, N 
and O, can penetrate each other by diffusion. Excitabil- 
ity would then be strictly restricted to a boundary region 
along the interface and this is what we basically see in the 
experiment. Using the diffusion values which have been 
used for quantitative simulation of the chemical wave pat- 
terns in Rh(110)/NO+H 2 we can estimate the diffusion 
length I at T=740 K for r=10 s with I = V2Dt resulting 
in I = 8 /im for N and 13/im for O [13]. The inset in 
Fig. la shows a dark boundary region of a few /im width 
which is consistent with the high WF of 1.1 eV of the 
c(2x4)-20,N phase [15]. 

For modeling the observed behavior we set up a di- 
mensionless 3-variable model for bistable/excitable me- 
dia which in 2D reads 

d t u =u — u 3 — v — S(u—u s )q 2 + d u Au + d uq Aq, (la) 

d t v =e(u + (3 - v) + d v Av, (lb) 
8 t q =(1 - q)(q - a)(q + 1) + 7 (1 - q 2 )(u - u s ) 

+ d uq Au + Aq, (lc) 

with diffusion constants, d u ,d uq ,d v > 0, parameters 
/3,7, 8 e R, e > and — 1 < a < 1. In short, using 
U = (u, v, q) with obvious notations we write 

d t U = f(U) + DAU. (2) 

The system (1) is composed of an excitable u, v- 
subsystem (FHN like) and a bistable g-subsystem (Allcn- 
Kahn or Nagumo equation), which has front solutions. 

The basic idea is that (i) through the interaction with 
the g-variable the u, w-subsystem is excitable only in the 
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vicinity of the front position (where q w 0), and that 
(ii) these localized excitations of the it, u-subsystem then 
push or pull the q- front. Since on surfaces the diffu- 
sion of the different species is not independent of each 
other, we include cross-diffusional terms which have to 
be symmetric according to Onsagers reciprocity relation. 
On surfaces cross diffusion arises (i) due to the vacant 
site requirement for diffusional hops and (ii) due to ener- 
getic interactions between coadsorbed species [16, 17]. In 
particular, the strong repulsive interaction between coad- 
sorbed oxygen and nitrogen shows up in a downward shift 
in the N 2 desorption maximum by about 100 K [18]. As 
will be shown below cross-diffusion becomes important 
for the nucleation of excitation pulses. 

Thus, we choose parameters j3,e in such a way that 
in the absence of q, i.e., for q = 0, the (u,v) ODE sub- 
system dt(u,v) = (fi(u, v,0),f2(u, v,0)) is excitable. Its 
unique ODE fixed-point (u s ,v s ) is given by u s — — /3 1 / 3 , 
v s = u s + /?. This fixed point is asymptotically stable 
and globally attracting, but for small e > rather small 
perturbations may lead to large excursions. 

For it = u s , or equivalently 7 = d uq = 0, (lc) is a 
standard bistable equation 

d t q = g{q) + Aq, g{q) = (1 - q)(q - a){q + 1), (3) 

i.e., the kinetics d t q = g(q) has the two stable fixed 
points u = ±1 and the unstable fixed point q = a. It is 
well known, that (lc) has travelling front solutions, e.g., 
q(x, y, t) = qf(x — cot), <?/(£) — > ±1 as £ — > ±00, in fact 
explicitly given by c = V2a and g/(£) = tanh(£/v / 2)- 
For a < (a > 0) fronts travel left (right), meaning that 
the +1 phase invades the —1 phase (resp. vice versa). 

Since the Laplacian is isotropic any orientation of 
fronts is allowed. As a consequence, (3) also has 
(smooth) V-shaped fronts qy, propagating with speed 
Ci = CQy/l + 1/m 2 , see Fig. 2 and [19]. 
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FIG. 2: Heuristics for V-shaped fronts of (3). 

Now considering the coupling between (la,lb) and (lc) 
in more detail we note that |d ug Ag| becomes large near 
corners of the front, and vanishes away from the front; 
thus (it, v) excitations originate near corners. On the 
other hand, the term — S(u — u s )q 2 in (la) makes the 
(it, v) kinetics less excitable away from the front, and thus 
excitations in the PDE (1) stay near the front. Finally, 
the term 7(1 — q 2 )(u — u s ) in (lc) has the effect that 



the excitations push or pull the g-front, as seen in the 
experiment. 

System (1) was integrated numerically in a domain 
f2 = [— L, L} 2 for various parameters using different initial 
conditions (IC) (it, v, <?)|t=o = ( u o, v 0i Qo) and boundary 
conditions (BC) . For the IC we are led by the experiment 
to consider "wedges" in q, e.g. 

*(*>») = {l x>x -m\y\> < 4 > 

where ±m € R are the slopes of the sides and xoGR 
represents the position of the tip. For (it, v) we choose the 
fixed point (u n , v ) = (u s ,v s ). Given an IC of the form (4), 
it is natural to integrate (2) in a moving frame £ = x — rjt 
with r] fts Ci(m) to keep the tip of the wedge away from 
boundaries, i.e., to integrate 

t U = f(U)+DAU + V d^U. (5) 

For the BC the problem then is that while planar fronts 
can be easily simulated with Neumann BC, for V-shaped 
fronts influences of boundaries on the fronts are difficult 
to avoid. Here we choose Dirichlet BC for (5), namely 

(u,v)\ an = (u s ,v s ) and q = ±1 on £ = ±L, 
q(Z,±L)=q f (Z-ti ). 

The latter fixes the front shape and position at the top 
and bottom boundary. 

For the IC and BC chosen above, we obtain the sim- 
ulation results displayed in Fig. 3. Excitations nucleate 
at the tip of the wedge and then travel along the front, 
pushing it back and forth. The chosen 7 = —0.05 < 
means that u > u s (it < u s ) pushes q down (up), such 
that here the excitations push back the frontline. The fir- 
ing process at the tip repeats for some time (essentially 
depending on the size of the computational domain) , and 
the process is accompanied by some overall reshaping of 
the wedge. Aside from boundary effects, this reshaping 
is determined by the following factors. The q- front does 
not fully recover its former position after a (u, v) pulse 
has passed. The tip of the wedge, where pulses nucleate, 
drifts to the right. To counteract this effect we chose 
77 = 3(7/4 (instead of 17 — c\ which without coupling to 
the (u,v) system would give a stationary tip position). 
As a consequence of decreasing |?7|, the unperturbed sides 
of the wedge drift to the left. The overall balance gives an 
almost stationary average front position up to t = 500. 
For t > 500 excitations that have emanated from the tip 
are reflected by the boundary. 

The behaviour in Fig. 3 is quite robust with respect 
to most parameters and ICs, including the opening angle 
of the wedge. A decisive parameter for the evolution 
is 7. For 7 = —0.2 the excitations push the front too 
strongly thus destroying the wedge by creating a bubble. 
For 7 = 0.1 the excitations pull the front too strongly 
thus flattening the wedge, see Fig. 4. 
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FIG. 3: Numerical integration of (1) in frame moving with 
speed rj — 3ci/4 = —0.15. Parameters d u = 0.09, d v = 
0.01, d uq = 0.1,/3 = 0.2,6 = 0.5, e = 0.03,7 = -0.05, a = 
—0.1. IC for q is the wedge (4) with xo = L/4, m = 1, ICs 
for (u,v) are (u s ,v s ). BC according to (6) with £o = — 3L/4. 
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FIG. 4: Same parameters and IC as in Fig. 3 except for 7. 



In summary, we observed excitability in a catalytic 
surface reaction which remained strictly localized at the 
interface of two domains of different adsorbates. The 
excitations were traveling along the interface in a pulse- 
like way, causing lateral displacements of the interface 
position. Mechanistically, the localized excitability can 
be traced back to the diffusive mixing of the two sep- 
arate adsorbates at the interface causing the formation 
of a mixed coadsorbate phase which is required to make 
the system excitable. The experimentally observed be- 
havior could be reproduced with a general dimcnsionless 



3-variable model which couples the excitability of a sub- 
system to the position of a frontline. The nucleation 
of excitations at corners of the front was explained with 
cross-diffusional effects which are very sensitive to the lo- 
cal front geometry. Similar dynamical behavior should be 
expected in all systems which (i) are essentially bistable 
in the sense that there are two asymptotically stable 
phases, but where (ii) diffusive mixing at the interface 
can locally change the dynamical behavior from bistable 
to excitable. 
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